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Abstract 

In supersymmetric models, scalar top quarks, or stops, generically have parity- violating 
couplings to top quarks. As a result, tops produced in stop decays should be polarized. 
In this paper, we will argue that this effect may be observable at the LHC with realistic 
integrated luminosities, provided that one of the stops is copiously produced and can 
decay to a top and a neutralino. Wc define the "effective" stop mixing angle, which 
determines the degree of top polarization, and discuss the prospects for a measurement 
of this angle at the LHC. If some information about the neutralino mixing matrix 
is available, this measurement can be used to constrain the mixing angle in the stop 
sector, one of the most important ingredients in assessing the naturalness of electroweak 
symmetry breaking in the MSSM. 



1 Introduction 



Supersymmetric (SUSY) extensions of the Standard Model (SM) are one of the most com- 
peUing candidates for physics at the electroweak scale [1,2]. Among many new particles 
predicted by these models, scalar partners of the top quark, the stops ii and t2, play a 
special role. The reason is that the top Yukawa coupling is by far the strongest coupling 
of the SM Higgs. The dominant quantum contributions to the Higgs potential arise from 
loops of the top quark and its superpartners. These contributions play an important role 
in the physics of electroweak symmetry breaking (EWSB) in this framework. For example, 
in no-scale unified models such as mSUGRA, they are responsible for triggering EWSB by 
driving the Higgs mass^ term negative, while all other scalar masses remain real. More gen- 
erally, the stop and top loops give a sizeable contribution to the parameters determining the 
EWSB scale and the Higgs mass spectrum. Thus, if TeV-scale SUSY is realized in nature 
and discovered at the LHC, determining the properties of the stop quarks will be crucial in 
understanding the mechanism of EWSB. 

In the minimal supersymmetric standard model (MSSM), sizable quantum corrections 
to the physical mass of the lightest CP-even Higgs are required to satisfy the LEP-2 lower 
bound on this mass [3] . The tension between this requirement and naturalness of the EWSB 
scale is known as the "LEP paradox" or "little hierarchy problem" [4]. Minimizing this 
tension requires a certain choice of the stop soft masses and the trilinear soft A-term [5], 
pointing to a "golden region" in the MSSM parameter space [6]. This region is characterized 
by fairly light stops (typically 300-400 GeV for ii and 500-800 GeV for ^2), and a large 
mixing angle 6t (defined as the rotation angle from the gauge eigenbasis to the mass 

eigenbasis ii, . The shape of the golden region depends only weakly on the non-top sector 
MSSM parameters (within broad ranges of those parameters), and within the unrestricted 
MSSM this hypothesis is compatible with other interesting applications of SUSY such as 
dark matter [7]. If SUSY is discovered at the LHC, it would be very interesting to test 
whether the golden region hypothesis is realized. 

There is a large body of literature on superpartner mass measurements at the LHC. 
(For example, Ref. [8] studied stop mass determination in a framework similar to the one 
considered here.) The goal of this paper is to propose observables that are directly sensitive 
to the mixing angle 6t- The basic idea is that this angle enters the stop-top- neutralino 
couplings. In particular, the amount of parity violation in these couplings is determined by 
9t (along with the neutralino mixing matrix). If the decay t — > tx^ is kinematically allowed, 
this parity violation results in non-zero polarization of the tops produced in this decay. As 
is well known, tops decay before they hadronize, so that the top polarization is reflected in 
the angular distributions of the top decay products [9]. Observing these distributions can 
therefore provide information on Ot- In this paper, we will show how such measurements can 
be performed at the LHC. 

The rest of the paper is organized as follows. In Section [2| we discuss the connection 
between the stop mixing angle and the polarization of tops from t —>■ tx^ decays in more 
detail, define the effective mixing angle which determines the degree of top polarization, and 



1 



and 



propose two observables which can be used to measure this effective angle. In Section |3| we 
perform a parton-level analysis of the feasibility of the proposed measurements at the LHC, 
for a particular benchmark point (chosen to be compatible with the golden region) in the 
MSSM parameter space. The analysis includes the leading SM backgrounds, and we discuss 
a set of selection cuts to suppress these backgrounds to manageable levels. We project the 
sensitivity of the LHC experiments to the effective mixing angle with 10 fb~^ of analyzed 
data. In Section |4} we briefly discuss some of the issues involved in inferring the actual stop 
mixing angle 9t from the measured effective angle. Finally, we present our conclusions and 
outline some questions for future study in Section [5] 

2 Polarized Tops from Stop Decays 

The stop mass terms in the MSSM Lagrangian have the form 

C = {tlXR} M^tL^tRf , (1) 

where 

^2 ^( Ml + ml + Au V^mt sin /3( A - /i cot (3) \ 

\ V^mt sin l3{At - /x cot /?) M^ + m'f + Au ) ' ^ ' 

1 2 \ 2 

- — -siv?6w\ cos 2/9 m|, = - sin^ 6*14/ cos 2/3 m| . (3) 

The mass eigenstates ti and ^2 are superpositions of the gauge eigenstates, ti and iji: 

ii = cos 9t II + sin 9t Ir , 

h = - sin Ot h + cos 6t in , (4) 
where the mixing angle is given in terms of the Lagrangian parameters by 

2y/2mtsm P {At - fi cot (3) 

Since parity is violated in weak interactions, weak couplings of the stops depend on the 
angle Of. Of particular interest to us is the stop-top-neutralino coupling, since the parity 
asymmetry in this coupling can lead to non-vanishing polarization of top quarks produced 
in stop decays. Since the top quark decays before it hadronizes, this polarization is, at least 
in principle, observable. The relevant vertex has the form 

g%Ux', {cose%P, + sme%Pn)t, (6) 

where j = 1 . . . 4 labels the neutralino mass eigenstates x^, related to the gauge eigenstates 
N = {B,W^H'„H'J by 

4 

X^^=J2N,kNk. (7) 

k=l 
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The effective mixing angles are given by 



tan^^^^ = 



V2 (f iV,2 + f iV,i) cos + VtN^^ sin Ot ' 
VtN^i sm9t + ^g'Nji cos Ot 



V2 ( |iV,-2 + |iV,i) sin Ot - VtNji cos Ot 

wfiere yt — V^mt/{v sin P). The main idea of this paper is that we may be able to get an 
unambiguous and fairly precise experimental measurement of one or more of the angles at 
the LHC, by measuring the polarization of top quarks produced in the decay t — > x'^t. If the 
neutralino mixing matrix is at least partly known from other measurements, this information 
can be used to extract (or at least constrain) Of. This information can in turn be used, 
together with the stop eigenmass measurements, to determine the stop-sector lagrangian 
parameters. 

At the LHC, stops can be directly pair- produced by strong interactions, in the processes 

PP ^ iii*. (9) 

Direct production of same-sign stop pairs is negligible. In addition, there may be a sizable 
sample of stops produced indirectly, namely in decays of other superpartners, particularly 
the gluino via g — > ti. (A brief discussion of the possibility of top polarization measurements 
in the gluino sample appeared in Ref. [10].) Those events lead to more complicated final 
state topologies in the detector, and vetoing such topologies can be used to separate the 
"direct" and "indirect" stop samples. We will focus on the direct stop sample in this paper, 
assuming that the contamination from the indirect sample, if present, is negligible. This 
has the advantage of simpler events and more robust predictions, since the rate and event 
topologies in the indirect sample depend on many more MSSM parameters. Some of the 
analysis techniques described here could be applied to the indirect sample as well. 

Once produced, stops will promptly decay. Possible two-body decay modes include t —>■ 
iX°) i ~^ ^X^) t — > W'^b, and t — > We are interested in the t — > tx° mode, which must 

be kinematically allowed and have a sizeable branching ratio for our analysis to apply. This 
decay is followed by t — > W~^b, and the W-boson then decays either hadronically (about 70% 
of events) or leptonically (about 10% for each lepton flavor). Angular distributions of the 
top decay products are sensitive to top polarization. For example, the angular distribution 
of the b quarks in the top rest frame has the form 

oc ( + 2\ {E^ + sin 20esm^) + - 2 ) cos 29^^ cos 9^, , (10) 



d cos 6*6 J 

where 9}, is the angle between the momenta of the b quark and the neutralino coming from the 
same stop decay as the top, and and are the energy and momentum of this neutralino. 
(See Appendix A.) In the case of hadronic W decay, the top rest frame can be reconstructed. 
The neutralino direction is unknown, since there is at least two missing energy particles in 
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each event. However, since stops are heavy, in the reaction (|9]) they are produced close to 
threshold, so that the stops in the direct sample tend to not have large velocities in the lab 
frame. If the stop is at rest in the lab frame, the top and neutralino from the same stop 
decay are back-to-back in this frame. Then, the angle Of, is the same as 6b, the angle between 
the b quark momentum in the top rest frame and the direction of the boost from this frame 
to the lab frame. The boost direction is simply the direction of the top momentum in the 
lab frame. The angle Of, is experimentally measurable. As we show below, the distribution of 
events in this angle provides a useful polarization analyzer, even though the effect is washed 
out somewhat by the motion of the stops in the lab frame. 

Another observable, useful for semileptonic top decays, is the angular distribution of the 
charged lepton. In the top rest frame it has the form 

oc + sm26esm^ + cos26cs cos 9i. (11) 

d cos 6i 

In semileptonic top decays, the missing neutrino does not allow for a reconstruction of the 
top rest frame. (Since there are other unobserved particles in the event, the two neutralinos, 
the neutrino energy and momentum cannot be reconstructed by imposing the W and top 
mass constraints.) However, one can define the "approximate" top rest frame in which 
Pb + Pi = 0) the angular distribution of charged leptons in that frame can be used as a 
polarization analyzer, as will be shown below. Once again, we will use di, the angle between 
the lepton momentum in the (approximate) top rest frame and the direction of the boost 
from this frame to the lab frame, as a stand-in for 6i. This preserves most of the polarization 
asymmetry due to the low lab-frame velocities of directly produced stops. 

Note that t and t decays can be combined in the polarization analysis, since the polar- 
ization effects are the same for the two samples as a consequence of CP conservation. 

3 Observability of Top Polarization at the LHC: a Bench- 
mark Study 

In this section, we demonstrate that top polarization in the stop pair-production sample can 
be observed and measured at the LHC. As with most effects in the MSSM, observability 
of top polarization, and the details of the analysis required to extract it from the data, 
depend on the choice of the MSSM parameters. We will perform a detailed analysis at a 
specific benchmark point in the MSSM parameter space, compatible with the golden region 
hypothesis of Ref . [6] . 
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Figure 1: Leading-order cross section for pp — > tt* at the LHC We set the factorization 
and renormahzation scales to n = mi and used the CTEQ611 parton distribution function 
set [15]. 

3.1 Benchmark Point and Signature 

We choose the following benchmark point (BP) in the MSSM parameter space (all parameters 
are specified at the weak scale): 

Ml = 100 GeV, Ma = M3 = 1 TeV; /x = 400 GeV, tan/3 = 10; 

"^(<5i,2) = "^(^ta) = "^('^1,2,3) = "^(-^1,2,3) = "^(6^2,3) = 1 TeV. (12) 

All phases, flavor-violating couplings and masses, and all A-terms except for At, are set to 
zero. The three stop sector lagrangian patameters, m((53), m{u'^, and At, determine the 
two stop eigenmasses, M(ti) and M(t2), and the mixing angle 9t- In our analysis, we dial 
the lagrangian parameters to keep the stop masses fixed at 

M(ti) = 340 GeV, M{t2) = 800 GeV, (13) 

while 6t is varied to illustrate the analyzing power of the top polarization measurement. 
Note that the chosen stop eigenmasses, tan j3 and /i are compatible with the MSSM Golden 
Region, for values of 6t sufficiently close to 7r/4 (maximal mixing). 

The cross section of direct stop-antistop pair production at the LHC, as a function of the 
stop mass, is shown in Fig. [!} At the BP, the cross sections ar^ 

(t{pp iiil) = 3.23 pb , 

a{pp ^ ^2^1) = 33 fb . (14) 

^All signal cross sections cited in this paper are evaluated at leading order, using CTEQ611 pdf set and 
setting the factorization and renormalization scale to the mass of the produced stop state. For the the ti 
background channels we set the factorization and renormalization scale to the top pole mass and for the 
VK+jets events we use /i^ = + X]jGts("^tr)' where rnjj. = + Pt^ summed over heavy quarks and light 
jets. We use the MadGRAPH/MadEVENT software package [11] to compute cross sections and generate events, 
except for the VF+jets background channels, where we use ALPGEN [12]. 
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Figure 2: One of the diagrams contributing to the direct stop production in proton-proton 



colhsion, followed by the decay chain leading to the final state (15). 



The direct sample is completely dominated by ti pairs, and we will ignore the small t2 con- 
tamination in our analysis. Moreover, all indirect stop production processes, such as gluino 
decays, are strongly suppressed at our BP; we will ignore those processes. The produced 
ii decays, with essentially 100% probability, via ii txi', all other two-body decays are 
kinematically forbidden. Thus, each event contains a top, an anti-top, and two neutralinos 
in the final state. Both hadronic and semileptonic top decays provide observables which 
can be used to infer top polarization. We will focus on the two angular distributions dis- 
cussed in Section [2} 9b in hadronic top decays and 9i for semileptonic tops. A priori, it is 
not clear which observable would provide better sensitivity in a realistic experiment: While 
the hadronic channel suffers from the finite resolution of the jet energy measurement and 
combinatoric backgrounds, the leptonic channel is limited by the inability to reconstruct the 
top rest frame in each event. We will explore this question by analyzing events where one 
of the tops decays hadronically and the other one semileptonicalljj^ allowing us to study 6^ 
and 6i distributions in the same sample. Thus, our final state is 



(15) 



where i = e, fi. An example of a Feynman diagram contributing to this final state is shown 
in Fig. |2} Other final states (two hadronic or two semileptonic tops) can also be considered, 
and we expect that some of the techniques discussed here will be applicable to those channels 
as well. However a detailed analysis is outside the scope of this paper. 



^For a recent study of the LHC reach for stops and spin-1/2 top partners in this channel, see Ref. [13]. 
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6.5 
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Table 1: Signal and Background cross sections (in pb), before and after cuts. Also listed are 
the total number of events simulated for our study (including the b-tag efficiencies), and the 
software package used to generate the events. 



3.2 SM Backgrounds and Cuts 

Before presenting the polarization analysis, let us discuss the SM processes contributing to 
this final state, and propose a set of cuts that strongly suppress their contributions without 
materially affecting the polarization asymmetries in the signal. 



Several SM processes contribute to the final state (15). The most important processes 
are listed in Table [T] For concreteness, we will give cross sections with ^ = The largest 
contribution by far comes from the tt channel, where one of the tops decays hadronically and 
the other one semileptonically. The semileptonic decays may produce /i" directly or in a r" 
decay; we list the two contributions separately, to emphasize that the cut efficiencies for the 
two contributions differ due to their different kinematics. Another important background 
process is 4j + W, where the invariant mass of three of the jets accidentally coincides with 
the top mass, and the W decays leptonically. In contrast, another obvious background, tiZ, 
where the Z decays invisibly, has a very small cross section (only about 0.7 pb) and thus we 
will not include it in the analysis. 

We generated and analyzed Monte Carlo event samples for the signal and each SM back- 
ground roughly corresponding to 10 fb"^ of integrated luminosity at the LHC. The simula- 
tions were performed using MadGraph/MadEvent (MG/ME) software package [11] except in 
the case of VT+jets background, where we used the ALPGEN package [12]. The MG/ME simu- 
lations included integration over the phase space of the fully decayed final state (8-body final 
state in the case of signal, 6-body for the tt backgrounds). In the matrix element calculations, 
small-width approximation was not used, but only the Feynman diagrams corresponding to 
the desired decay chain, such as the one shown in Fig. [2| were included. All polarization 
effects in the decay chains of interest are fully taken into account in this approximation. To 
verify this, we compared the simulated distributions in the angles 9h and 9i for the signal to 



■^For the signal and the ti background channels, the cross sections with i — e^, e+, /i^, /i+ arc identical. 
For the VF+jets background channel, the cross sections with a negatively charged lepton are suppressed by 
about a factor of 2 compared to the ones shown in Table [T] 
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Figure 3: Comparison of the analytic predictions and events generated by 
MadGraph/MadEvent for the angular distribution of the b (left) and the charged lepton (right) 
in the top rest frame. Here 6b {9i) is the angle between the momenta of the h quark (charged 
lepton) and the neutralino coming from the same stop decay as the top. The solid (red) his- 
togram corresponds to the Monte Carlo events with cos26'efr = — 1 and the dashed (red) line 
corresponds to the analytic prediction. The dashed and dash-dotted (black) lines correspond 
to the MC distribution and analytic prediction, respectively, for cos2^eff = +1- 



the analytic predictions in Eqs. (10), (11). The results are in excellent agreement, as shown 



in Fig. [3] Our analysis is restricted to the parton level; effects of hadronization, showering 
and intial state radiation are not included. To roughly model the detector response to jets 
and electrons, we introduce a Gaussian smearing of their energies according to [14] 



50% 



© 3% 



GcV 



10% 



0.7%. 



(16) 



GcV 



We do not apply smearing to muon energies, since the effect is small. To model the acceptance 
of the detector, we apply the following acceptance cuts at the generator level: 



• pT,j > 25 GeV, rj < 4.0 for each jet; 

• Pt > 20 GeV and rj < 2.5 for the charged lepton; 

• AR{ji,jk) > 0.4 for each jet pair, and Ai?(jj, i) > 0.4 for the charged lepton and each 
jet (z,fc = 1...4). 

The results of our analysis are summarized in Tables [T] and [2j The first column of Table [T] 
lists, for reference, the total cross section of each process. For the W+ jets processes, the 
listed cross sections are computed imposing a cut ptj > 20 GeV on light (non-bottom) 
quarks and gluons. The second column of Table [T] lists the cross sections to produce the final 
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state (15) in each channel, including the above acceptance cuts. At this point, the signal 
is completely swamped by the tt background, with S/B 0.01. To extract the signal, we 
impose the following selection criteria: 



• At least one of the jets must be 6-tagged. We assume energy-independent 6-tag prob- 
abilities of 50% for a true &-jet, 10% for c-jet, and 1% for light quark and gluon jets. 

• Large missing transverse energy: ^ > 125 GeV. 

• cos (j){p^, ^) < 0.7, where 0(p^,^) is the opening angle between the transverse mo- 
mentum of the charged lepton and the missing transverse momentum. This cut elimi- 
nates events where all missing energy comes from the leptonic decay of a highly boosted 
W. 

• Hadronic Top Reconstruction (HTR): min |m3j — m^'''^'^| < 5 GeV, where rrisj is the 
invariant mass of three jets, and the minimum is over the four possible triads. In our 
study, the true top mass mf^^ is taken to coincide with the value used in the MC 
generator, 172.5 GeV. The jet triad for which the minimum is achieved is identified 
with the hadronic top. This cut suppresses iy-|-jets backgrounds. 

• Semileptonic Top Veto (STV): Assuming that all missing energy comes from a massless 
neutrino, we reconstruct the candidate neutrino four-momentum p'^ using the trans- 
verse measurement and the conditions (p^)^ = 0, {p^ +p"Y — where is the 
momentum of the jet which does not belong to the triad identified as the hadronic top. 
There are generically two solutions, Pi2- We then form sj'^ = {p^ +P'iYi and demand 
miuj \ \/sf- — mw\ > 40 GeV. This cut suppresses the tt backgrounds, since the re- 
jected kinematics corresponds to missing energy coming exclusively from semileptonic 
top decays. 

• Separation Cuts (SC): AR{i^,X) < 1.5 , where AR{i^,X) is defined as the AR 
between the charged lepton and the jet closest to it in A_R; and 0.8 < A0(^, X) < 1.3, 
where A(j){^, X) is defined as the azimuthal angle between the pr vector and the 
object closest to it in 0. The signal and background distributions motivating these 
cuts are presented in Figure |4j 

The efficiencies of these cuts for the signal and each of the background samples are collected 
in Table |2| and the cross sections after cuts for each process are listed in the third column 
of Table [T} With these cuts, we obtain 

S/B = 2.5; S/Vb = 36 {dnt = 10 fb"^). (17) 

Thus, we conclude that the stop signal is easily observable above the SM backgrounds after 
the selection cuts are imposed. We can then use this signal-dominated event sample for the 
polarization analysis. 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 

AR x) (MET, x) 

Figure 4: Normalized distributions in AR{i^,X) and A0(^,X), in the signal and back- 
ground samples satisfying all other acceptance and selection cuts. (See the text for the 
definition of the plotted observables.) The solid (red) line corresponds to the signal and the 
dotted (gray) line corresponds to the combined backgrounds. 
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38 
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37 


24 
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Table 2: Selection cut efficiencies, in %. (Each cut is applied to the sample satisfying the 
acceptance cuts and all selection cuts listed above it in the text.) 



Before proceeding, let us comment on the reliability of the Monte Carlo predictions for 
backgrounds, such as the ti simulation used in this study. Convincing interpretation of 
any excess over the SM background as contribution from new physics, in a situation where 
S/B ^ 1, would require precise understanding of the assumptions used in the background 
predictions. Since the LHC will produce a very large sample of ti pairs, the MC generators 
used to simulate this background can be precisely calibrated with data. This calibration can 
be done using samples which are not expected to suffer any SUSY contamination, such as 
the events with two reconstructed hadronic tops and small It can then be applied to 
samples where SUSY can contribute, e.g. events with large ^ and a charged lepton which 
are our main focus here. It would be interesting to understand quantitatively the expected 
accuracy of such a validation procedure; such an analysis is however outside of the scope of 
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Figure 5: Angular distributions of events in the angle 6h. The different contributions cor- 
respond to (from top to bottom): signal (yellow), 4j + W~ (black), 2j + 2b + W~ (white), 
ti{fi~) (gray), tt(r~ /i~) (light red). The event numbers correspond to 10 fb~^ integrated 
luminosity at the LHC 



this paper. 



3.3 Polarization Analysis 

In Section[2| we identified two observables, the angles Ob and 6i, which are potentially sensitive 
to the top polarization in the direct stop sample. The angle 6b is the angle between the b jet in 
the jet triad identified as the hadronic top, and the direction of the hadronic top momentum, 
in the rest frame of the hadronic top. To define 6i, we first define the approximate rest frame 
(ARF) of the semileptonic top, by adding the three-momenta of the charged lepton and the 
jet that does not belong to the reconstructed hadronic top. Then, 9i is the angle, measured 
in the ARF, between the charged lepton and the direction of the ARF's motion with respect 
to the lab frame. 

To measure the angle db, we need to find which of the three jets in the hadronic top is 
the b jet. One way to do it would be to simply demand that one of the jets be 6-tagged. 
An alternative method, which only relies on kinematics, is to compute the invariant masses 
of the three jet pairs inside the hadronic top, and to find the pair whose invariant mass is 
closest to mw The jet that does not belong to that pair can then be identified as the b jet. 
We found that in our sample the kinematic method has a higher efficiency than the simple 
6-tag, and thus we will use this technique to compute 9b in each event. 

In Figures |5] and [6| we show the distribution of events in our Monte Carlo sample (signal 
and background, passing the above acceptance and selection cuts) in cos 6';, and cos 9i. The 
event numbers correspond to 10 fb~^ integrated luminosity at the LHC. Each figure shows the 
distributions for three values of 9es: from left to right, 9es = 0,7r/4, 7r/2. The asymmetry 
due to top polarization is clearly visible in the distribution. One should note that, even 
though we were careful to choose selection cuts that have as little effect on the hadronic 
asymmetry as possible, the cuts do have a slight preference for the events with cos^;, < 0, 
resulting in small asymmetries for the backgrounds and the signal at the parity-conserving 
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Figure 6: Angular distributions of events in the angle 6i. The different contributions cor- 
respond to (from top to bottom): signal (yellow), 4j + W~ (black), 2j + 2b + W~ (white), 
ti{fi~) (gray), ti{r^ (light red). The event numbers correspond to 10 fb~^ integrated 

luminosity at the LHC 



point, ^eff = 7r/4. Likewise, in the leptonic asymmetry, our use of the ARF instead of the true 
top rest frame leads to non-zero asymmetry in these cases. The important point, however, is 
that even though the measured asymmetry is not identical to the "primordial" asymmetry 
due to top polarization, there is still a simple one-to-one map between the two. To quantify 
the effect, we define two forward-backward asymmetries (FBAs): The hadronic FBA is given 
by 

[!l-j\)d cos 9, 

while the leptonic FBA is defined as 



du 
^ dcosO), 



{lo-J-i) d COS 01 



da 



A% = ^ . (19) 

' d cos 01 



The statistical uncertainty of the FBA is given by 

^Afb = 2 ; 20 

where N'^ = d cos 9i j^ff and is the number of events in the opposite hemisphere. 
The combined asymmetry is defined as 



/I combined _ / /ihad /ll*^PN ( Ahad A^'^P \ (0T\ 

-^FB — \^FB ~ ^FB) ~ \^FB ~ ^FB) ) \'^^) 

where the second term signifies the average over the cos26'eff values. For our MC event 
sample, we obtain 
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Figure 7: Leptonic, hadronic, and combined forward-backward asymmetries, as a function 
of the angle O^s- The error bars indicate statistical errors for 10 fb~^ integrated luminosity. 
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where the errors are statistical only, corresponding to 10 fb^^ worth of data at the LHC We 
conclude that, at the level of our analysis, polarization effects in the stop direct sample are 
easily observable. 

If the masses of ti and Xi are measured independently, the measurement of the FBAs 
can be used to determine 6'cfT- We studied the dependence of the hadronic and leptonic 
FBAs on ^eff by generating MC samples of the signal for 10 values of 6es, keeping all other 
MSSM parameters fixed at their BP values. The results are shown in Fig. [7} The error 
bars indicate statistical errors for 10 fb~^ integrated luminosity. As expected, we observe an 
approximately linear relationship between the asymmetries and cos2^eff: 

±0.12 cos20efr - 0.41 , 
-0.037 cos2^eff- 0.047, 

-0.16 cos2^eff- (22) 

These formulas can be used to estimate the precision of ^efr determination. Note, how- 
ever, that the coefficients in these formulas depend on the stop and neutralino masses, and 
would need to be recalculated if these parameters differ from the BP values. The errors 
in the experimentally determined masses would need to be taken into account in the 9es 
determination. 

Among the effects not included in our parton-level analysis, arguably the most important 
one is the absence of extra jets emitted as initial state radiation (ISR). Such jets would 
introduce a combinatoric background, making it more difficult to identify hadronic tops and 



.lep 
^FB 
A had 
^FB 
A combined 
^FB 
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Figure 8: Leptonic and hadronic forward-backward asymmetries for p^j > 40 GeV, as a 
function of the angle ^gff- The error bars indicate statistical errors for 10 fb~^ integrated 
luminosity. 



potentially washing out the forward-backward asymmetries. The pr distribution of the ISR 
jets is sharply peaked at zero, since most of these jets are soft or coUinear with the incoming 
beams. An additional lower cut on the jet pt (above the acceptance cut of pxj > 20 GeV 
assumed in the above analysis) may be necessary to suppress such jets. To estimate the 
sensitivity of our results to such a cut, we have repeated the analysis requiring ptj > 40 
GeV, keeping all other cuts the same. The results are presented in Fig. |8| It is clear that the 
hadronic asymmetry becomes essentially useless in this case, while the leptonic asymmetry is 
only marginally affected. The main effect of the additional ptj cut on the leptonic asymmetry 
is in fact simply due to the reduced statistics, which can be overcome in time by increased 
integrated luminosity. Thus, we are optimistic that the leptonic asymmetry would survive as 
a useful observable even after the ISR and other showering, fragmentation and detector effects 
are included. Of course, a full detector-level analysis is needed to confirm this conclusion. 



4 From Top Polarization to Stop Mixing. 

If the effective angle 9es is measured using the technique described in the previous sections, 
the actual stop mixing angle can be found by inverting Eq. ([s]): 

ytNi4 cos ^eff - V2 f f iVi2 + ^Nn) sin O.s 

tan Ot = ^ . (23) 

ytNii sin ^efi + ^g'Nu cos ^efr 

Note, however, that several issues arise in this conversion. First, the top polarization ob- 
servables measure cos2^eff, which introduces a two-fold degeneracy in the determination of 
Ot- Secondly, in some cases, the mapping from cos2^efT to 6t introduces large errors even if 
cos 29cs is precisely known: for example, this can occur for pure-bino neutralino if cos 2^eff is 
close to one, as seen in Fig. |8j Finally, to determine 6t, one needs to know the composition 
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Figure 9: cos20cff vs. the stop mixing angle 9 for three different neutrahno mixing matrices. 
The sohd (black) curve is our benchmark scenario, the dotted (gray) curve corresponds to 
the case where the lightest neutralino is pure bino and the dashed (red) curve corresponds 
to pure higgsino. 



of the lightest neutralino in terms of the gauge eigenstates. For example, in the cases of pure 
bino and pure higgsino, drastically different values of 9t would be inferred from the same ^eff 
(see Fig. [s]). In the case when the neutralino is predominantly a bino, it is especially im- 
portant to independently constrain the subdominant contributions, since the mt and WHt 
coupling constants are larger than the Btt coupling, so that even a small higgsino or wino 
admixture can have a significant effect of O^q. Constraining the neutralino mixings at the 
level required to make the 6'efr — > Ot conversion possible may in fact require substantially 
more integrated luminosity than the top polarization measurement itself. In spite of all 
these complications, this strategy seems to offer our best hope for the direct determination 
of the angle 9t in the unconstrained MSSM framework, and is worth pursuing if SUSY is 
discovered. 

5 Discussion and Conclusions 

In this paper, we argued that top and anti-top quarks produced in stop decays in the MSSM 
are generically expected to be polarized, and that observing and measuring this polarization 
can provide information about the stop mixing angle. We proposed observables sensitive to 
top polarization, and illustrated the potential of the LHC to observe this effect, by performing 
a detailed study for a single benchmark point in the MSSM parameter space. We conclude 
that for the chosen benchmark point, the polarization effect should be observable with 10 
fb~^ of data, and a useful measurement of the "effective" stop mixing angle 6^^ can be 
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obtained. We also briefly discussed how the measurement of ^es can be used to constrain 
the actual stop mixing angle 6*4, one of the crucial parameters in assessing the naturalness of 
the EWSB in the MSSM. 

The most important limitation of our analysis is that it is performed at the parton 
level, and only a very rough modeling of finite detector resolution by smearing jet and 
electron momenta was included. It is very important to perform a detector-level analysis 
including the effects of showering, hadronization and fragmentation, as well as detector 
effects. In particular, as we already remarked, ISR jets could contribute serious combinatoric 
backgrounds not included in this study, particularly for low-p^ jets. These backgrounds can 
be suppressed by tightening the minimum ptj cut, and we have verified that at least one of 
our observables (the lepton asymmetry) appears to be rather insensitive to such tightening 
(see Figs. [7| [s]). Still, it is important to analyze these issues quantitatively with a detector- 
level analysis. 

Acknowledgments 

We are grateful to Johannes Heinonen, Jay Hubisz, Mihoko Nojiri, Matthew Reece, 
Christian Spethmann and Gilad Perez for useful discussions. MP would like to thank the 
Institute for Physics and Mathematics of the Universe (IPMU) in Tokyo, and MP and AW 
would like to thank the Kavli Institute for Theoretical Physics (KITP) in Santa Barbara, 
where parts of this work were completed. AW is grateful to the Aspen Center of Physics 
where parts of this work were done. This research is supported by the NSF grant PHY- 
0355005. 

Note Added: The main idea of this paper and some early results have been presented at 
the KITP Conference "Anticipating Physics at the LHC" in June 2008 [16]. Related ideas 
have also recently appeared in Ref. [17]. 

A Observables Sensitive to Top Polarization 

The top quark decays before it hadronizes, and the information about its helicity is reflected 
in the distributions of the decay products. Consider the process ti x^t, followed by the 
decay of the top quark t W^b. The matrix element for this process is proportional to 

u{p^) (cos 6esPL + sinOesPR) y^^Ui{pt)ui{pt)l^PLu{pb)tt,{pw) ■ (24) 

i=± 

It is straightforward to calculate this matrix element in the coordinate system where the top 
quark is at rest, and the neutralino momentum is along the +z axis. The non-zero matrix 
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elements in the helicity basis are 

M{x+:Wo) = f COS 9es y/ Ey^ - + sin O^s V^x+ Px) ' — ? : 
A< (X+, VTj.) = (^cos 6'eflf ^E^-p^ + sin O^s \/^x + Px) ■ cos - , 

COS 



M{x~,Wq) = (cos^effA/E^T+P^ + sin^effA/KT^p^ 
M (X- , W^±) = (cos e^s y/E^+p^ + sin ^eff y/E^-p^j • ^2 sin - , 



(25) 



where E^ and p^ is the energy and momentum of the neutralino, 9 is the angle between the 
b quark momentum and the z axis. The differential cross section is proportional to 



+ 2^ {E^ + sin26'effm^) + - 2^ p^ cos26'eff cos6' . 



(26) 



The presence of a term linear in cos 6 indicates parity violation, and measuring the size of 
this term provides information about the effective mixing angle ^eff- 

An analogous calculation for the leptonically decaying top results in a matrix element 
squared proportional to 

\Mf (X E^ + sin 29esm^ + p^ cos 26'efr cos §1 . (27) 

Note that in both cases the magnitude of the neutralino momentum in the top rest frame is 
fixed because the neutralino is the final state of an effective two body decay. 
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